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Abstract 

We theoretically analyze diffusion trajectories of an anisotropic object moving on a two dimen¬ 
sional space in the absence of an external field. In determining diffusion parameters associated with 
the shape anisotropy, we devise a measure based on the gyration tensor, and obtain its analytic 
expression exactly. Its efficiency and statistical convergence are examined in comparison with the 
fourth cumulant of particle displacement. We find that the estimation of diffusion constants based 
on the gyration measure is more efficient than analysis adopting the fourth cumulant. 
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I. INTRODUCTION 


Particles suspended in fluid move along random trajectories. Study on so called the Brow¬ 
nian motion has provided fundamental understanding of macroscopic diffusion phenomena 
led by collective dynamics of microscopic constituents. Theoretical framework to describe 


rm 


the Brownian motion was proposed by Einstein [l| , which is experimentally verified by Per 


2[. The Brownian trajectories reflect medium properties, shape and size of particles. 


Hence, analyzing diffusive motions is a useful way to examine the physical properties of 
solution and colloid suspension, and its applications extend to the realm of chemistry and 
biology, let alone physics 

Although the original theory by Einstein describes translational motion in terms of dif¬ 
fusion constants, theory was extended to rotational motion of an ellipsoid js, 9|, and an 
object with shape anisotropy has acquired considerable attention over decades for both in 


free space 


lJl5| 


and under an external potential held 


16| . There exist abundant exam¬ 


ples of anisotropic objects such as proteins and microtubles in nature, man-made colloidal 
particles [l7|, ll8|, and nanotubes 19l-l2l|. It is worthwhile to mention a few diffusion char¬ 


acteristics of anisotropic particles. Since object more easily diffuses in the direction of its 
major principal axis, trajectory of an anisotropic particle becomes elongated along the pre- 


Q- 


ferred direction [HJ. This asymmetric trajectories in turn allow us infer the shape of a 
particle, but it can only be observable within time interval while orientation of the particle 
is preserved. Also, the statistics of diffusion trajectory is not Gaussian having higher order 


cumulants which contain the information associated with shape 13l-ll5|. 


It is only recent that these features could be experimentally characterized 
of art experiment was enabled by single particle tracking technique 


141. The state 


22|, and demonstrates 


a potential possibility to acquire shape information from diffusion trajectories. Yet, the 
practical application is limited by a few factors. First is to realize high temporal resolution. 
If time resolution is longer than orientation correlation time, diffusion trajectory becomes 
indiscernible from that of spherical or point particles in the laboratory frame, and infor¬ 
mation related to shape anisotropy is totally erased. This issue would be partially solved 
by recent developments to enhance the resolution less than orientational correlation time of 


membrane proteins of cells, which is typically order of 10 ~ lOO/x sec 23|, [2^. More crucial 


is the fact that it is experimentally infeasible to identify initial angles of nanoscale objects. 
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and one has to consider averaged qnantity over an ensemble of randomly oriented parti¬ 
cles. In conventionally used second cumulant of displacement, terms associated with shape 
anisotropy vanish upon averaging over randomly distributed initial angles. We thus needs 
to consider higher orders as relevant quantities to the trajectory analysis in relation to the 


shape factor 


13l-ll5l|. Since statistical convergences of higher order cumulants are usually 


poor, their statistical behaviors and errors are still to be examined. Furthermore, natu¬ 
ral question arises: What is the best higher order quantity to pursue efficient and reliable 
estimation of diffusion constants and shape anisotropy? 

The purpose of this study is to propose a new measure, based on the gyration tensor, 
and suggest a systematic way to read diffusion constants of anisotropic objects from particle 
position trajectories. Considering motions of anisotropic objects moving on a two dimen¬ 
sional space in the absence of external force, we analyze the fourth order fluctuation of 
particle displacement both analytically and numerically by Langevin dynamics simulation. 
In particular, we obtain the analytic expression of the gyration measure, which completely 
agrees with Langevin dynamics simulation results. We compare the statistical convergence 
and hnite sampling error of the gyration measure with previously suggested fourth moments 
and hnd that analysis adopting gyration tensor is most efficient. This paper is organized as 
follows. In Sec. II we introduce our system and its governing equations of motion, where 
we illustrate an example of the diffusive trajectory, and discuss the properties of the second 
cumulant of particle position. The higher order quantities are examined in Sec. Ill, where 
we introduce the fourth order quantity obtained from the gyration tensor and obtain its 
analytic expression exactly. In order to test the efficiency of those measures, we evaluate 
the standard errors of the fourth cumulants and the gyration tensor based measure and 
explicate their errors in the estimation of diffusion constants in Sec.IV. Finally in Sec. V, 
we discuss a possible scheme to implement least error measure and the properties of other 
fourth order measures. Summary will follow in Sec.VI. 

II. SYSTEM 


We consider a Brownian particle of arbitrary shape, moving on a two dimensional space. 
In the absence of external force, the translational motion of the Brownian particle can be 
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FIG. 1: Two dimensional diffusive motion of a highly anisotropic object with D±^ = 0 (the upper 

panel) in comparison with a spherical particle (the lower panel). Trajectories are obtained by 

Langevin dynamic simulation, where total time elapse is with time resolution (At) chosen 

to be equal to the angular relaxation time 1/Dq for the whole trajectory. The local trajectories 

(see the blow ups) are acquired for finer time resolution At = Q.QI/Dq, and display distinctive 
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features depending on object shape. 









described by the positional Langevin equation as 

n (t) = f dt'ii {t'), i = x,y, ( 1 ) 

where (r^yTy) denotes the center of mass coordinate, and the noise components and 
have zero means and statistical properties given as 

= (*-«')■ ( 2 ) 

The noise amplitude is determined by diffusion matrix Dij{6). Here we dehne the particle 
orientation coordinate 6 to be an angle between x axis and one of eigenvectors of diffusion 
matrix Dij, and ■)e denotes the noise average for a hxed orientation 9. On the other 
hand, the equation of rotational motion of the particle is similarly given by 

e{t)= f dt'ie it'), (3) 

Jo 

where is white noise with = 0 and = 2 D 0 S{t — t') with Dq being the 

rotational diffusion constant. Here we dehne the particle orientation coordinate 9 to be an 
angle between x axis and one of eigenvectors of diffusion matrix Dij. Note that Dij is a 
symmetric matrix which is diagonalizable to have eigenvalues, say, and D± (diffusion 
constants along major and minor axis, respectively), and in laboratory coordinates, can be 
represented as: 

Dij (9) = (l/2)[H(5ij + S{az cos29 + axsm29)] (4) 

with A = {D\\ + D±) and S = {D\\ — D±), and being the Pauli matrices. Three diffusion 
constants. A, S and Dg, fully characterize motion of the object, where nonzero S indicates 
anisotropic motion, resulting from shape anisotropy. In order to verify the analytic results to 
be derived shortly, and in order to explicitly evaluate hnite sampling errors, we employ the 
Langevin dynamics simulation in which position and orientations are updated by integrat¬ 
ing randomly generated noises, Eqs ([I]) and ([3]), satisfying respective huctuation dissipation 
theorem. Time is discretized in the simulation with unit time step, say At. Length of At 
should be much smaller than characteristic time scale of the phenomenon dealt in the simu¬ 
lation, which is the orientational relaxation time, 1/Dg, for this case. We set At = O.Ol/D^ 
for generating trajectories presented in Figure. 1. Positions of the object are recorded at 
every time At, which is time resolution of the observation. 
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Figure 1 displays typical translational trajectories of a particle described by Eqs. ([1]) and 
(j^. where we compare the diffusion trajectory of an ellipsoid to that of a sphere. The 
trajectories are obtained by choosing a time resolution equal to the orientation relaxation 
time, At = I/Dq, and the trajectories in the marked areas are acquired by higher time reso¬ 
lution, At = 0.01/Dq. While the overall shape of diffusion trajectories of a sherical particle 
remains self similar upon changing the time resolution, the trajectory of the anisotropic 
particle traced with the shorter time resolution appears different from the coarsely resolved 
trajectory, signaling the correlation between angular and translational motion. 

In characterizing stochastic trajectories, it is instructive to examine the moments of a 
particle position at a given time. In the system of our interest, the hrst moment is identically 
zero due to the noise property. Using Eqs. ([I]) and ([2]), one hnds the second moment, 

dti [ dt2{^x{ti)^x{t2))e 
Jo 

= At cos 29 (1 - , (5) 

where the average is performed over ensemble with a hxed initial angle 9. The second 
moment of y coordinate of the particle trajectories can be obtained by replacing S with —S. 
For an isotropic object (S' = 0), the second moment describes well the normal diffusion and 
the position dispersion is perfectly symmetric, as indicated by the aspect ratio, = 

1. For S' 7^ 0, the diffusive behavior in the long time limit [Det 1) is similar to the normal 
diffusion with diffusion constant A/2, and the aspect ratio of trajectories is almost unity. If 
observation time is short compared to the angular relaxation time (Det 1), the diffusion 
follows the behavior, {r‘/. y)e ~ (A ± Scos29)t, dependent on the anisotropy factor S and 
the initial angle, and leads to the aspect ratio deviated from unity. 



III. FOURTH ORDER MEASURES 


In order to read S and Dq from the second moment given by Eq. ([5]), one needs to 
realize an ensemble of particles at a hxed initial angle, which can possibly be done if one 
detects the initial angles of particles, and collect a number of trajectories of particles which 
have the identical initial orientation. However, it is difficult to resolve the orientations 
of microscopic objects in single particle tracking experiment 1^ so that experimentally 
available are quantities averaged over randomly distributed initial angles. Meanwhile as 
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is apparent from Eq. ([5]) the average of the second moment over initial angle 6 nullihes 
the second term containing factors S and Dq. Therefore, it is necessary in most practical 
situations to examine higher moments in order to obtain the physical constants of objects 


from the initial angle averaged trajectories. To the end, the authors of Ref. 
the fourth cumulant of particle positions as 


14| obtained 


(C4 (i)> = YeS 

with C^i^t) = — 3 (r^)^. Here the average is also performed over initial angles, and the 

diffusion constants S and Dq appear as relevant parameters to the fourth cumulant. 

There certainly exist alternative ways to measure fluctuations in quartic order of the 
stochastic variable. For example, one can consider the fourth moment of translocation, (r^). 
Following the calculation scheme in Appendix A, we hnd, 

(r^) = SAV + (1 - ,-».<) . (7) 

L>0 ZUg 

Subtracting the hrst term given from the second moment of r leads to a measure, M 4 = 

4.4_2(r2)2: ^ 

(M.) = + e-“>“ - 1 ). ( 8 ) 

For short time measurements {Dgt 1), (M 4 ) increase quadratically in time and in the 
long time measurement {Dgt ^ 1), it converges to linear function of time. Such asymptotic 

behaviors are determined by S and Dg. This is also the case for {C 4 ) because of a relation 


(^ 4 ) = 3(M4)/8. 

Another fourth moment measure is gyration tensor, which is useful in identifying the 


spatial distribution of particles 
positions: 


25 


27|. The gyration tensor is constructed based on particle 


Gn — 


2 m 


N 

*,i=i 


^ijUij Uij 


( 9 ) 


where N is the total number of particles. Here distance between the Ah and the jth particle 
are denoted as Xij = — Vxj along the x-axis, and yij = Vy^i — Vyj along the y-a.xis. The 

gyration tensor, Eq. ([9]), has two non-negative eigenvalues, say Ai and A 2 which act as 
shape identihers and do not depend on the choice of coordinate. Addition of eigenvalues, 
A+ = Ai -|- A 2 gives the squared radius of gyration which corresponds to the size of particle 
distribution. The difference of the eigenvalues, A_ = |Ai — A 2 I, signihes the anisotropy of 
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FIG. 2: Quartic order measures for analyzing diffusive motion of an anisotropic object {S/A = 1) as 
a function of observation time t in units of the orientational relaxation time where numerical 
results (points) from Langevin dynamic simulations are averaged quantities acquired from 10® 
simulations. In comparison, analytic expressions (lines), Eqs. ([6]) and (llip . are plotted, showing 
an excellent agreement with numerical results. Straight lines are the asymptotic behaviors given 
in Eq. (fT^ . 


the distribution; A_/A+ vanishes for an isotropic distribution, and becomes close to unity if 
the distribution is highly anisotropic. 

For the present problem to analyze particle trajectories, we consider that and Vy^i in 
Eq. on]) represent the position of a particle at discrete time ti, that is, {t/) = r^^i, and 
ry{ti) = j. For the fact that A_ quantihes an anisotropy of distributions of Tq, j’s, one may 
take A^ as a quartic order measure. However, the average of A^ contains terms independent 
of S and Dq, which should be subtracted to dehne a relevant measure. We introduce a 


gyration measure as Mg = A^ — /a with /a given by 



( 10 ) 





where (A+) = At/3. After some lengthy but straightforward calculations, we obtain exact 
expression for Mg (see the details given in Appendix B) : 


(Mg) 

F{x) 


F{2Det) 

1920A)2 

15e“^'^(l + xY + 15(x^ — 1) — + 4:X^, 


( 11 ) 


which is again function of S and Dg, and acts as another relevant measure for analyzing the 
diffusive motions of anisotropic objects. 

The fourth order measures discussed so far entail factors related to the shape anisotropy 
in their respective forms. In particular asymptotic behaviors, which can be summarized as 


{X)D^g/S^ 


a{Dgt)‘^ , Dgt <C 1 
(ttDgt , Dgt S> 1 , 


( 12 ) 


allow us to extract the diffusion constants of anisotropic object. Here the coefficients for X = 
C 4 , M 4 , Mg are a = 3/2,4,1/9 and /3 = 3/4, 2,1/15, respectively. This long time and short 
time behaviors are displayed in Fig. 2, where we compare numerical results obtained from 
Langevin dynamics simulation with the analytic results. It can be seen that the two methods 
yield consistent results. This consistency could be obtained only for sufficient repetitions of 
simulation (10® simulations were performed in producing Fig. 2). If the number of repeated 
simulations (or experiments) is insufficient, the estimate of X would be randomly deviated 
from the analytic results. Therefore, in practical estimation of the diffusion constants from 
a finite number of measurements, the statistical convergence of the measures and also finite 
sampling errors in estimating the diffusion constants come into question. 

Before proceeding, let us make a remark on the effect of a marker position. In single 
particle tracking experiments diffusion trajectory of an object is obtained by tracing a point¬ 
like marker allocated on the object. In our consideration, the marker is assumed to be placed 
at the center of mobility. Even if the marker position is displaced from the center by £, 
one can show by repeating the similar calculation that the long time behavior of Eq. (IT^ 
remains unaffected, while only the short time behavior is modified: The effect of the off- 
center distance on the short time behavior can be completely described by replacing S with 
S — Dgi'^ in Eq. flT^ . From a dimensional analysis, one thus finds that this correction would 
be negligible when when Dj_) ^ with L the size of the elongated object. 
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IV. STANDARD ERRORS AND FINITE SAMPLE BIAS 


In order to investigate the statistical convergence, we consider the variance of the mea¬ 
sures, 

AA = ^{X^) - (A )2 . (13) 

Although it is the very basic idea of statistics, we briefly mention the statistical meaning of 
AA. Suppose that one has a data set of Aj’s with Aj obtained from the ith measurement 
(or simulation), and the total number of data, n, is relatively large so that one can apply 
the central limit theorem. Then, the bias of the hnite sampling average 
from the true average (A) can be written as 

A,, - (A) ~ XXj^ . (14) 

From the above equation we hnd that the number of samplings required to reach the desired 
accuracy, e* = (A„ — (A))/(A), amounts to 

AA/(e*(A)) . (15) 

Therefore, AA in relative to the average (A) quantihes the hnite sampling error and also 
required number of measurements. 

In Fig. [3l we present the relative huctuation, AA/(A), as a function of observation time. 
We observe that it increases linearly in time, and has hnite ohset at t = 0, and hence 

AA/(A) ~ k {Dot) + C. (16) 

Here the slop k and the ohset C depend on the choice of measure and the shape anisotropic 
factor S/A, as presented in the inset of Fig. 3. We also perform the detailed analysis of 
statistical errors of the gyration measure which indeed conhrms Eq. flT 6 |) and allows us to 
analytically evaluate k and C. Over the entire range of observation time, AA/(A) for 
A = 04 is bigger than that for the gyration measure A = Mg. Moreover, O 4 increases 
more rapidly than Mg, and its statistical convergence in the long time limit should be much 
poorer than Mg. This observation together with Eqs. fHTD and f[T5]) implies that if one 
choose to use O 4 rather than the gyration measure, estimation of S and Dq according to 
the limiting behaviors Eq. (IT^ from hnite number of experiments can be more errorneous, 
and requires more samplings. In this respect it appears that more efficient analysis can be 
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FIG. 3: Relative fluctuations, Eq. m, of the quartic order measures for various S/A values as 
a function of Dgt, where the average is taken over 10® simulations. Relative fluctuation linearly 
increases in time, and it has an offset at Dgt = 0, showing the behavior given in Eg. (1161) . The 
inset shows the dependence of the offset and the slope on the shape factor S/A and the choice 
of measnre. Symbols represents values obtained from simulation while analytic expressions for C 
for Mg (solid line), k for Mg (dashed line), C for (^4 (dotted line), and k for C 4 (dash-dot line) are 
plotted for comparison. The fluctuations for C 4 are more significant than Mg for a given S/A, and 
for a given measure, it is more pronounced for smaller S/A. 

done by the gyration measure. However, because of AX/{X) linearly increasing in time, in 
using the gyration measure it is desirable to choose the observation time window, t < 
where tm is large for the long time behavior to set in, and yet not too large to bring in 
unwanted fluctuations. It should also be noted that the relative fluctuation becomes more 
significant if we consider less anisotropic particles (smaller S/A). Note also in Appendix C 
that the magnitude of k and C are enhanced as S/A decrease. This is because for small S 
the measure itself becomes small, while {X"^) remains hnite. Therefore, it is expected that 
a large number of experiments should be repeated in order to probe weak anisotropy. 
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Let us now demonstrate more explicitly how the discussed behaviors of the variance are 
reflected in the estimation of diffusion constants. We evaluate hnite sampling errors in S and 
Do, and present the results in Fig. 4, where and represent values extracted from 
hnite sample average of X of sample size n, X^-, with X = and X = Mg. In obtaining 
SD)^ we measure the coefficient of of using the short time behavior in Eq. flT^ . In 
order to ensure Dg <^1 condition, the slope is evaluated at t = with unit time step 

Ar = \Q~^/Dg . In the upper panel of Fig. 4, we plot the bias, — S'!/S' with S being 
a true value, as a function of sample size. As can be seen, the hnite sample bias decreases 
with n, and follows well l/y/n behavior presented by the straight line. For the smallest 
sample size {n = 10 ^), the error is roughly order of ten percent, and yet C 4 yields the error 
a few times larger than Mg, which is due to its large huctuation as previously stated. The 
statistical huctuations of C 4 , more signihcant in the long time limit, cause serious problem in 
the estimation of Dg. As shown in the lower panel of Fig. 4, the bias ranges roughly from 1 
to 10 , and moreover, data points do not converge well along the line showing l/y/n behavior 
due to its poor convergence. On the other hand, for the gyration measure, although the 
error is large, it is order of magnitude less than the case of O 4 , and more importantly, it 
can further be reduced by increasing the sample size, as suggested in the bias well saturated 
into I/a/u. 


V. DISCUSSION 


In the preceding sections, we show that in analyzing trajectories of an anisotropic particle 
the gyration measure is more efficient than the fourth cumulant of displacement. Before 
concluding our work, it may be instructive to discuss a formal procedure to implement the 
least error measure and explain involved complications. One may write the generalized 4th 
order measure as. 


A 


dtirai,iC, 


Q;iQ;2Q:3Q;4 


(t, ti, t2, ts, t^) , 


i=l 


( 17 ) 


where represent the position of particle along ai = x,y direction at time t*. The mea¬ 

sure with least error should be invariant under permutation of x and y as well as interchange 
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FIG. 4: Finite sampling errors as a function of sample size, where and are the estimated 
values from measuring the asymptotic behaviors of The upper panel shows the bias of 
from its true value S. Both C 4 and Mg result in the bias decreasing with the sample size (l/y/n 
lines are added for guides to the eye), but the bias is larger for C 4 in the entire range of sample 
size. In determining D^\ we evaluate the slope of at t = The bias from 6*4 is order of 

magnitude larger than Mg, and convergence t(\'^/y/n behavior is not yet achieved for the sample 
sizes in the plot range. We also choose the observation time t = 2IDq, although not shown here. 

























of time variables, and therefore, there are only three possibilities for Ca^a 2 a 2 ,a 4 , to give 

‘^ = 17/ dti[Ci{xiX 2 X 2 ,x^ + yiy 2 y 2 ,yi) (18) 

i=i 

+ C2 {x^ysVi + yiX2X3X4) + C3 {XiX 2 y 3 y 4 )], 


where Ci = Cxxxx, C 2 = Cxyyy, and C 3 = Cxxyy Arguments of Ci are omitted and the sub¬ 
indices of X and y represent the indices of time variable. Note that X in general contains 
terms irrelevant to shape anisotropy. With these anisotropy irrelevant terms subtracted from 
A, we obtain a measure, say, X. Relative fluctuation determining the statistical error of a 
measure X is given by i?[{C'}] = XX/{X) with AX defined in Eq. (ir3|) . In order to acquire 
X having least error, we need to build an equation for a function set (( 71 ,( 72 , C 3 ) which 
minimizes the functional E'[{(7}]. Although stating the problem may seem straightforward, 
there exist several hurdles to overcome. One of those is without a function set of (7’s 
specified, actual calculations are hardly doable. Besides, even if we choose a special set of 
{(7}, extremum condition for E'[{(7}] leads to a complicated integral equation for (7 [2^. We 
therefore cannot resort to this generalized scheme mentioned above. 

Alternative way to find an efficient measure is to consider possible candidates of specific 
forms, and directly compare their statistical errors with one another, which is the strategy 
chosen in this work. In addition to ( 74 , M 4 and Mg, we take other measures such as 




dtix\, 


n = n 

i=l 

^ fi rt 
i=l ^ -'U 


dti 


dti 


{xjxl + yly^) , 

{X1X2X3X4 + 1/1I/22/3I/4), 


and examine their statistical properties. Here without calculation details, we only present ob¬ 
tained results that in the long time limit the errors of the measures listed above also linearly 
increase in time as XXi/{X/} = kiDyt with different coefficients, kiR^ = 8\/M/(3\/5) ~ 
7.35, k 2 R^ = 2\/l838/(3\/^) ~ 4.83, and = 80/(9\/3) ~ 5.13. These values are larger 
than kR/^ = 4\/^/(3\/7) ~ 2.81 for the gyration measure (see Appendix C), which sub¬ 
stantiates that the gyration measure outperforms not only the conventionally used fourth 
cumulant but also other measures considered above. 
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VI. SUMMARY 


We considered the thermal motion of anisotropic objects in two dimensions, which is 
determined by three diffusion constants -Dy, D^, and Dq. For the goal to determine the 
diffusion constants by analyzing diffusion trajectories, we investigated measures of quartic 
order in particle position such as the fourth cumulant ( 6 * 4 ), the translocation measure (M4), 
and the gyration measure {Mg). The translocation measure is independent of initial orien¬ 
tation of a particle, and related to the fourth cumulant as {M 4 ) = 8 (U 4 )/ 3 . The measure 
{Mg) is obtained from eigenvalues of the gyration tensor, which identify the shape of par¬ 
ticle trajectories. We exactly obtained analytic expressions for the quartic order measures, 
and confirmed that numerical simulations give consistent results with the exact solution. 
It was found that the asymptotic behaviors of the measures are determined by the value 
of S' = D\\ — Dx and Dq, as given in Eq. ([12]). We further evaluated the variance of the 
measures in order to test their statistical convergence. It was shown that the relative vari¬ 
ance increases in observation time and has hnite offset at short time limit. The offset and 
the slope of the increment is small for {Mg), in comparison with or M 4 , suggesting the 
efficiency of the gyration measure. We examined finite sampling error in the estimation of 
diffusion constants using and Mg, and demonstrated explicitly that Mg yielding smaller 
error outperforms the other measures. We also showed that even using Mg, estimation of 
Dq from the long time behavior can be severely biased due to the property of the variance 
increasing in measurement time, which is also the case for the evaluation of not only Dq but 
also S for a particle with weak anisotropy. 
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Appendix A: Evaluation of M4 


We consider a displacement of a particle from an origin, r = and evalnate its 

fourth moment, 

{r^) = {Tl) + 2{rlrl) + {rl) . (A.l) 

In this appendix, we let to represent an average over a random variable x which 

includes ^x,^y,^e and the initial angle. According to Eq. ([I]), the positions can be expressed 
as the integrated noises, and the terms in Eq. flA.ip can be written as 

4 .t 


) = 

\ ai a2' 




dtk{i ai OL2 01-2 (^ 4 ) ) 


with Ofc = x, y. Since the noises are Gaussian variables, one can apply Isserlis’ theorem 
(or Wick’s theorem) to write the average of quartic product of the nosies in terms of the 
averages of paired noises: 




(A.2) 


2)(e C1:2,3^Q:2)4) 2)(e ai , 3 ^ 0 : 2 , 4 ) 


where abbreviations = iai,k are used. Using the property of the translational noise 

Eq. ([2]), we obtain 

= Vij{ei,e^)5{ti-t2)6{U-tA) (A.3) 

= 4A.(^i)/^,,(^3) + 8A,(0i)A,(03) , 

where 6i = 6{ti). Summing up the results obtained so far, and using the diffusion matrix 
components Dij given in Eq. (jlj), we reach 

= [ dt' f + 4A^cos20(t',f")] (A.4) 

Jo Jo 

with Q(t',t") = 6{t') — 0{t"). We are now left with average over angular noise and initial 
angle. This angular average process should be performed on the cosine term in the above 
equation. For the Gaussian property of we have the following identity: 

^ (A.5) 

Plugging this relation for n = 2 into Eq. flA.4D . integrating with respect to time leads to the 
fourth moment of translocation, Eq. ()7|) . 
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Appendix B: Evaluation of 


Let us first evaluate the average of A+ which is the sum of eigenvalues of Gat in Eq. 

. N 


x^ = 


2 iV2 


[(^,2 G,j) + f'yj) ] • 


(B.l) 




We take ra,i to represent the position of particle along a = x,y direction at time U. The 
averages of the components appearing in Eq. fIB.lIl are given by Eq. (|5]) to lead 

N 


= (B- 2 ) 

In the continuum limit of time, we can replace the summation with an integration as 

N 

(B.3) 


N-^ V ^ f dt' , 

i=i 


and obtain 


lim (A+) = ^ [ dt' [ dt"A \t' — t"\ = -At. 

Ai^o' e_L /„ ' '3 


(B.4) 


The subtraction of the eigenvalues yields a lengthy expression of A?_ given in terms of 
combinations of quartic multiples of particle positions at different times. We give a short 
hand writing of the expression, 


N 


A^ = 




Ir*.**' + + (x <-> a)], 


(B.5) 


{*fc} = l 


where the summation denotes a multiple summation, q h P’ ^ indi¬ 

cates terms produced by permuting x and y of the terms given in the square braket. The 
summand, is given only in terms of x coordinate: 

7 ^(^) = t ■ t ■ t ■ t ■ t ■ t 

'x ' x,ll' x,l2 x,ll' X,12' X ,13 ^ ' x,ll' x,l2' X, 13 ' x,l4 1 

(‘2') ^ 2 ) 

and the mixing term Vx Vy contains factors quadratic in both x and y coordinate: 

(2) (2) o 2 I 2 2 

I X f y — x^i^ f y,i2 • y,i3 ' ’ I y,i2 • x,Z3 • yM ' x,ii ’ y,i2 

I Orv^ . rp .nr, . rp - Ann .nr, .nr, .nr, 

' x,ii' y^ii’ x,i 2 ' y,2 x^ii' y,ii' x,i 2 ' y^i3 • 

Analytic expression for gyration measure are obtained in the continuum time limit where 
the approximation in Eq. flB.3ll can be applied. We define (Aig) = hmAt^.o (A?_), and 
re-express it in concise form as 

(^ 9 ) = (n 7 / dtkraj^itk) ) , (B. 6 ) 


{«fe} 


\fc=i 
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where au = x,y. Summation combinations appearing in Eq. fIB.Sp . 

and the corresponding coefficient is given by C^a}- In obtaining (Aig), we first evaluate the 
average of quartic moment of position variables: 


n ’■«. (*t))=/ii / 


(B.7) 


^ '•ife 

\k=l / \fc=l“'° 

Upon using the Isserlis’ theorem, the average of quartic moment of noises can be written as 


) = (1.2)(3,4) + {1,3){2,4) + (1,4)(2,3) 


\k=l 


where {k,i) = Upon averaging over the translational noises, we have 




and the time integration of these quantities in Eq. flB.7p can be done, for which we define 

/•min(4,q) 


4>k,e = 2 




Here min is smaller value between tk and ti. Then, the average of quartic moment of 

spatial variables is given by 

= ( 01 , 203,4 + 01 , 302,4 + 01 , 402 , 3 )^^’^° ) 

\k=l / 


where the average over the angular noise can be derived through a relation for the Guassian 
property of ^g\ 

in/S.6(t)^irn/S.6{t')\^Q _ — 

O / o 


with A6{t) = 9{t) — Oq. This gives the average of the fourth moment of position at different 
times in Eq. flB.7j) for a given set of a^’s. Then Aig in Eq. fIB.bp can be evaluated by 
performing the integrations and summing over all the sets of a^’s corresponding to the 
terms in Eq. fIB.Sp . This process is lengthy but straightforward to reach 

+ 60{Det)^ - 80{Detf + 128{Detf]. (B.8) 


Defining Mg to eliminate the shape irrelevant term, that is. Mg = Aig — f a with = 
4(Ht)^/45, we obtain Eq. fITT]) . 
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Appendix C: Error coefficients of the gyration measure 


The statistical error of the gyration measure is evaluated, which leads to explicit expres¬ 
sion for k and C in the long and short time regime, respectively. In calculating the variance, 
AX for the gyration measure X = Mg, we use an identity, 

AM, = - (JV4)2 = - (M,y . 

Given (Mg) by Eq. flB. 8 D . we are left with evaluating (Mg). In long time regime, we only 
keep terms of highest order in time t and obtain 

(M^) = 128A^(3L^ - I 2 LIL 2 + SLl + 8 L 1 L 3 - 8 L 4 

-\- 2 L ‘^02 — 4T2O2 -l- O2 -\- 2O4) -|- 

Here symbols Li and Oi are defined below through an auxiliary variable, pj = dt'C,{t') 

with C,{t) being a Gaussian stochastic variable which satisfies {((t')) = 0 and (C(^)C(^O) = 

Li = 

L 2 = 

L 3 = 

L 4 = 

Oi = 

O 2 = 

O 3 = 

O 4 = 


n 


2 = 1 
3 


Ylliipm) {P2P3) = 


2=1 

4 


15 


T—r 17 

(P1P2) {P2P3) (psPi)= 


2=1 

5 


315 


(P1P2) {P2P3) (PsPi) (PAPb) 


2=1 


62 

2^ 


h (piPi) = -t 
2 

Yl^ipm) {P 2 P 1 ) 


= -L 


2 = 1 
3 


Ylhipm) (P2P3) (p 3 Pi) 


2 = 1 
4 


= —r 
15 


17 


(P1P2) (P2P3) (PsPa) (PiPl) = 


2 = 1 


(G.l) 
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where an integral operation is defined as Jj = dti/t. One can find that the relative variance 
or the statistical error is given by 


AMg/{Mg) ^ 


4 


(C.2) 


which confirms the behavior of Eq. flTbll linearly increasing in time and correspondingly gives 
the slope fc, 


, 4 1 

“ 3 V 


(C.3) 


with R = S/A. 

Next, let us evaluate the offset of AX/{X), C, in short time limit. Considering that 
orientation is almost fixed during short time interval, we find after some tedious algebra 
that (A4^) is expanded as: 


(M^) ~ 128A\MI)a + 128A^S\MI)as + IGS^M^s (C.4) 


with the coefficients given by 


(Ml)A = 3 Lf - 12Ly2 + 8 Ll + 8L1L3 - 8L4 

+ 2LIO2 - 4L2O2 + OI + 2O4 

AS = 9E3 — 30 L^L 2 + 8L2 + 28E1E3 — I6E4 

- 6L?Oi + I6L1L2O1 - I2L3O1 + LlOl 

- 2L2OI + LIO2 - 2L1O1O2 + OIO2 

- 4 L 1 O 3 + 4 O 1 O 3 + 4 O 4 

{Mfjs = 9 Lf - 24 LIL 2 + l6Ll - 12 LlOi + I6L1L2O1 
+ WLjOj - 8L2OI - ALiOf + Of + 12 Lf 02 

- I6L2O2 - 8L1O1O2 + 40 f 02 + 40 ^ . 


Inserting equations in (IC.ip into Eq. flC.4p . we obtain 

608 , 2608 


(-M ) = 


14175 


44 ^ /l2c2.4 I o4.4 

+ 14175 * +25* 


and further using that (Mg) ^ S'^t^/9 in the short time limit, as given in Eq. flTSp . we finally 
reach 


C = hmAM„/(M„) 

2 /l24 + 582/?2 + 98i?4 x 

“ sV 7 y’ 


(C.5) 
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The analytic expressions of k and C, respectively given in Eq. (IC.3P and Eq. fIC.Sj) . are 
consistent with the simulation result shown in inset of Fig. 3. 
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